ZÝ Aug 2009 


arXiv:0901.3818 


Residual effects of heavy sparticles in the bottom quark Yukawa coupling: 
a comparative study for the MSSM and NMSSM 


Wenyu Wang!, Zhaohua Xiong!, Jin Min Yang? 


1 Institute of Theoretical Physics, College of Applied Science, 
Beijing University of Technology, Beijing 100020, China 
2 Key Laboratory of Frontiers in Theoretical Physics, Institute of Theoretical Physics, 
Chinese Academy of Sciences, Beijing 100190, China 


If the sparticles are relatively heavy (a few TeV) while the Higgs sector is not so heavy (ma is 
not so large), the Higgs boson Yukawa couplings can harbor sizable quantum effects of sparticles 
and these large residual effects may play a special role in probing supersymmetry at foreseeable 
colliders. In this work, focusing on the supersymmetric QCD effects in the hbb coupling (h is the 
lightest CP-even Higgs boson), we give a comparative study for the two popular supersymmetric 
models: the MSSM and NMSSM. While for both models the supersymmetric QCD can leave over 
large residual quantum effects in hbb coupling, the NMSSM can allow for a much broader region of 
such effects. Since these residual effects can be over 20% for the hbb coupling (and thus over 40% 
for the ratio Br(h — bb)/Br(h — 7*77)), future measurements may unravel the effects of heavy 


sparticles or even distinguish the two models. 


PACS numbers: 14.80.Cp,12.60.Fr,11.30.Qc 


I. INTRODUCTION 


Supersymmetry is a prime candidate for new physics 
beyond the Standard Model (SM). Among various su- 
persymmetric models the most extensively studied is the 
minimal supersymmetric model (MSSM) fi]. Another 
popular supersymmetric model, which may be equally or 
more attractive compared with the MSSM, is the next- 
to-minimal supersymmetric model (NMSSM) | since it 
can solve the p-problem and alleviate the little hierar- 
chy. These models will soon be put to the test at the 
LHC. Thus, their phenomenological study is important 
and urgent. 

Although the most convincing evidence of supersym- 
metry is the detection of sparticle (sfermions, gauginos 
or Higgsinos) productions, the indirect probe through de- 
tecting the quantum effects of virtual sparticles in some 
measurable interactions will play a complementary role. 
If the sparticles are relatively heavy (say above a few 
TeV) and hence cannot be directly detected at the LHC, 
the indirect probe through quantum effects could be im- 
portant. For this end, the Higgs boson Yukawa interac- 
tions may play a special role since they can harbor siz- 
able quantum effects of heavy sparticles when the Higgs 
sector is not so heavy (i.e. ma, is not so large). The 
dominant quantum effects of sparticles are from the su- 
persymmetric QCD interaction and in the literature the 
calculations have been performed in the MSSM for such 
supersymmetric QCD effects in the Higgs boson Yukawa 
couplings a, Fi and the associated Higgs production pro- 
cesses at the LHC [5] The studies in the decoupling 
limit with heavy sparticles showed that for a light Higgs 
sector (ma is not large) the supersymmetric QCD can 


leave over large residual quantum effects in both the de- 
cay h — bb Z) and the productions at the LHC [5]. Given 
the popularity of the NMSSM, it is necessary to extend 
the study to the NMSSM. This is the aim of this work. 
In this work we will focus on the supersymmetric QCD 
effects in the hbb coupling (h is the lightest CP-even neu- 
tral Higgs boson) and perform a comparative study for 
the two popular supersymmetric models: the MSSM and 
NMSSM. Such a study is interesting for two points: 


(i) The study of the NMSSM can accommodate the 
study of the MSSM, and in some limit the NMSSM 
results can reduce to the MSSM results. One can 
envisage that the supersymmetric residual effects in 
the MSSM can be magnified in the NMSSM. So the 
future measurements of such effects may be useful 
in telling the difference of the two models. 


(ii) Compared with the MSSM, the NMSSM may pre- 
dict a different tree-level hbb coupling and different 
loop contribution. In the NMSSM the mass ma- 
trix and mixings of the Higgs bosons are enriched 
and thus the components of the lightest CP-even 
Higgs boson h are different from the MSSM val- 
ues. The residual supersymmetric QCD effects in 
the NMSSM may be larger and more interesting. 


Note that supersymmetry is a decoupling theory and all 
low-energy observables will recover their corresponding 
SM predictions when the mass scale of all supersymmet- 
ric particles (including the masses of sparticles and ma) 
take their heavy limits. The large residual quantum ef- 
fects of sparticles in the Higgs Yukawa couplings happen 
only in case that sparticles are heavy but m4 is light. If 


both sparticles and my, take their heavy limit, the resid- 
ual effects of supersymmetry do vanish. Since so far such 
a split scenario (with light Higgs bosons and relatively 
heavy sparticles) remains possible, we should check its 
phenomenological consequence. 


II. CALCULATIONS 


We start our analysis by recapitulating the basics of 
the NMSSM. In the NMSSM a singlet Higgs superfield S$ 
is introduced and the Higgs terms in the superpotential 
are given by 


Ss, (1) 


where A, and Ha are the Higgs doublet superfields, and 
A and « are the dimensionless constants. Note that there 
is no explicit u-term and an effective -parameter is gen- 
erated when the scalar component (S) of S develops a 
vev S: Heff = As. The corresponding soft SUSY break- 
ing terms are given by 


A, 
— AAS Ha: Hy — EA + h.c.. (2) 


So the scalar Higgs potential is given by 
Vr = |AHa- Hu — K9? |? + JASI? (Hal? + |Hul?) , (3) 


2 
Vo = & (\Ha?|Hul? — |Ha: Hal?) 


2 2 
+ 
+TP (|Ha? -|Hul?)” , (4) 
Vot = malHal? + mi |Hul? + m3|S/? 
= (A>àSHa -Hu + AnS? af hc.) (5) 


where gı and g> are the coupling constant of Uy (1) and 
SU,(2), respectively. So we can see that in the limit of 
vanishing A, k, A, and with the input of the effective 
Leff, the NMSSM can reduce to the MSSM. 

With the vevs vu, vq and s, the scalar fields are ex- 


panded as 
Va + Pa +t 
ae a (6) 
H} 
Hu = a Fal E , (7) 
S = stot+ié. (8) 


The mass eigenstates can be obtained by unitary rota- 
tions of interaction states, e.g., for the CP-even neutral 
mass eigenstates (h, Hı, H2) we have 


(h, Hy, Hz)” = V2U" (Qa, bu, o)”. (9) 


So the lightest CP-even neutral Higgs boson h is com- 
posed by 


h = UË V24 + UEV 2b, + UE V20. (10) 


In the MSSM, h can be decomposed in the same way 
except that without the singlet component, U¥ is a 2 x 2 
matrix and can be parameterized in terms of a mixing 
angle a. 

The sbottom squared-mass matrix in the NMSSM is 
same as in the MSSM with p replaced by perf. In the 
basis of (bz br), the squared-mass matrix is given by 


M? mX 
2 = L b#*b 
Mg Fe M? ) i (11) 
where 


Xp = Ap + perf tanG, 
M} = Mô +m + MZB — Qis) cos 26, 
MR = Mọ, + mi + MZQosiy cos 2/8. (12) 


with sw = sin fw, tan 8 = v,/va, IÈ and Qs being re- 
spectively the isospin and electric charge of the b-quark, 
Mo and Mz, being the soft breaking masses, and Ap 
being the soft breaking trilinear coupling. The sbottom 
mass eigenstates (bi, b2) are obtained by the unitary ro- 
tation of the interaction eigenstates: 


br B by _ ( cos® — sin 0; by 
(# )=20(# \=( 7 cos 0; bo (13) 


where the unitary matrix Zp parameterized by a mixing 
angle 6; diagonalizes the mass matrix in Eq. (11). 

The coupling of sbottoms with the singlet component 
o of the Higgs boson h comes from the F-term of the 
superpotential, given by 


g2™p 


— ~*~ )u, V2ob* br + h.c.. 14 
2mw cos 8 vuV2obj bn + É i 


In terms of the mass eigenstates of sbottoms, the vertex 
hb;bj takes the form 


g2™b 


Vie 5 vi 
where c’/ is only for the NMSSM (or equivalently set 
U{4 = 0 for the MSSM), while c? and ct! are present for 
both models. They are given by 
fem- m? cos? 3 Os 


ilj 
a Fei] 2828 


ij 
Cg 


2m? s?y cos? 8 


2 = 
j | ag 3Mp 


| 2328 
+4 (ZB ZH + ZBZB) , (16) 
ij m2, sin 28 2 soli 
2.2 ú 
Mz sy Sin 28 9; 595 
—— BZ 
3m; D/D 


thess (ZEZ + ZB ZB) , (17) 


V/2\my sin 8 
92 


ci = (Z5Z3 +2328). (18) 


Now we calculate the SUSY QCD corrections to the ver- 
tex hbb. At tree level it takes the same form in both the 
MSSM and NMSSM, given by 


r gomp UE 
V- = iA 19 
hob i 2mw cos 3 a) 
The one-loop SUSY QCD corrections come from the 
Feynman diagrams shown in Fig[I]where we do not show 
the self-energy loop of b. In our calculations we use the 


(a) (b) 


FIG. 1: Feynman diagrams for the SUSY QCD corrections to 
the vertex hbb. 


on-shell renormalization scheme [a and take the external 
b and b quarks on shell to get the effective vertex. With 
the corrections the effective vertex takes the form 


Vaos = Ving (1+ AY + A"), (20) 


where AY denotes the vertex correction from Fig[I{a) and 
A* is the counter term from the renormalization of mp 
and the wave functions of b and b. They are given by 


a a UE UA 

AY = (c+ SB 4 BC 
37 Ç tgga +ga? 
x [mo (ZZB + ZBZ) Cu + mg (2827 
+ZB Z ) Co] (mê, mè, m3, m,m? mē, ), (21) 

c mg. 

At = -e f Be ican [Bann 
—Bo(ms, m3, Mg, )] — 2m5 [B1 (mb; m3, mi, ) 


+Bi(ms,mj,mj, )] — 2memg sin(20z) 
x [Bo(ms, m5, m? ) - Bibs ma (22) 


where a, is the strong coupling constant, mg is the gluino 
mass, and Bo 1, Boa and Co,11 are the scalar loop func- 
tions [s] which can be calculated by using LoopTools [d]. 


Ill. NUMERICAL RESULTS 


In our calculation we use the package NMSSMTools 
for the mass spectrum and the rotation matrix of 

the Higgs fields (to get the corresponding MSSM results 
we take the limit of very small values for A, « and A,, 
and results are checked by using the package FeynHiggs 
{11}). To study the decoupling limit of SUSY parti- 
cles, we assume all the soft breaking mass parameters 


(M Mzp:MgAb) and the parameter fers are degener- 
ate, which are collectively denoted by Msusy. Then 
in the MSSM the SUSY parameters are (Msusy, Ma, 
tan 3), while in the NMSSM there are three additional 
parameters (A, k, A,). In our calculations we scan over 
these three additional parameters in the ranges 


~0.5<A,K < 0.5, — 500 GeV < A, < 500 GeV (23) 


Note here we did not consider a large A or «k. Theo- 
retically, the requirement of perturbativity up to some 
cut-off scale will set upper bounds on A and « at weak 
scale (if the cut-off scale is chosen to be the GUT scale, 
a stringent bound à? + k? < 0.5 is obtained (al). Phe- 
nomenologically, a large A or « will incur stringent con- 
straints from current experiments (see the last reference 
in i2). 

Before displaying the numerical results, we make some 
clarifications regarding to our numerical calculations: 


(1) In our calculations we used the package NMSSM- 
Tools which considered the loop corrections (es- 
pecially the stop/sbottom loops) to the effective po- 
tential, the masses and mixing angles of the Higgs 
bosons. These corrections are important and can- 
not be ignored. 


(2) Since we used the package NMSSMTools in which 
the Higgs mass matrices are diagonalized numer- 
ically, we did not use any approximate unitary 
transformation for our calculation. 


(3) For the b-quark mass mp, in the package NMSSM- 
Tools (and thus in our calculation), it is taken as 
the running mass mp(Q) (we take Q = Msusy), 
which means that the sizable QCD loop effects 
are taken into account. Actually, although both 
the tree-level coupling V,°.; and the one-loop SUSY 

QCD contributions ôV, are proportional to my 

and thus very sensitive to the value of mp, the rela- 

tive correction effects 5V;,,5/V,’,; (displayed in our 
numerical results) are not so sensitive to the value 
of Mp. 


Now we present some numerical results. To see the gen- 
eral feature of the corrections, we first switch off the ex- 
perimental constraints on the parameter space and in 
Figs[2] and B] display the resluts of the relative correction 
effects 


Asecp = 


OV pos 
p= AY + At. (24) 
hbb 
We see that for a light ma the SUSY QCD with a large 
Mgsusy can leave over sizable effects in the coupling hbb. 
But as m4 gets heavy, such residual effects of SUSY QCD 
become small in magnitude, showing the decoupling be- 
havior of supersymmetry. Compared with the MSSM 
results, the effects in the NMSSM can vary in a much 
broader region. While the corrections are always nega- 
tive in the MSSM, in the NMSSM the corrections can be 
both negative and positive. 
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FIG. 2: SUSY QCD corrections to hbb coupling versus ma in 
the MSSM and NMSSM for tan 8 = 20. The scatter plots are 
for the NMSSM while the curves for the MSSM. 


0.2 


m, = 100 GeV m, = 200 GeV 


m, = 500 GeV 


1000 2000 1000 2000 1000 
Mgysy (GeV) Mgysy (GeV) Mgysy (GeV) 


FIG. 3: Same as Fig. but versus Msusy for fixed values 
of ma. 


2000 


In order to figure out in what areas of the NMSSM pa- 
rameter space the one-loop SUSY QCD effects on the hbb 
coupling are sizable, we present a set of sample points in 
Table [I| We see that for the SUSY QCD loop effects on 
the hbb coupling to be sizable in the NMSSM, the ratio 
UJUE ( UR and Uff are respectively the components 
of d, and ġa in h, as defined in Eq[I0) plays the key role 
for the following reasons. As shown in Eq.(Q24), Asgcp 
is composed of two parts: the counter-term part A% and 
the vertex-loop part AY. We found that in most of the 
parameter space allowed by the LEP constraints, A‘ is 
negative and dominant in size, which is independent of 
how h is composed of. Whereas, AY is positive and can- 
cel A% to some extent. The size of AY can be enhanced 
by the ratio U#Z/U#{, as shown in Eq.(21). Although the 


TABLE I: Some sample points in the parameter space of the 
NMSSM, in which the SUSY QCD loop effects Asgcp can 
be of very different size. The mass parameters are in unit of 
GeV. Here we fixed Ma = 200 GeV and tan 8 = 20. With 
such fixed parameters and Msusy varying from 400 GeV to 
2 TeV, the value of mp in the MSSM varies from 106 GeV 
to 119 GeV while the corresponding values of Asgcp in the 
MSSM are about —11%. 


À k Akg Msusy UR/Uñ UR/UÄ mn Ea 
0 
0.106 -0.215 446 1788 16.2 -0.023 117.4 -4.1 
-0.116 0.180 333 1575 14.7 0.042 116.3 -6.0 
0.133 -0.256 393 1145 14.8 -0.046 114.7 -6.1 
0.123 -0.164 460 1568 14.1 -0.062 116.1 -6.8 
-0.103 -0.072 -294 1544 10.1 0.118 114.8 -12.0 
0.107 0.073 -175 1831 9.7 -0.100 115.2 -12.4 
0.126 0.101 -426 1738 7.3 -0.069 114.8 -15.5 
-0.106 -0.123 -377 1401 5.4 0.024 115.3 -18.0 
0.110 0.130 -324 1419 4.1 -0.017 114.9 -19.7 
0.106 0.406 -304 623 0.5 0.002 99.4 -24.5 


ratio UB JUE can also enhance the size of AY, its effect 
is suppressed by the smallness of UŽ and A in cd (in our 
scan we found that both \ and UË are small in order 
to satisfy the LEP constraints encoded in the NMSSM- 
Tools). Therefore, as U/U# gets large, AY becomes 
more sizable and, due to its cancellation effect, the total 
correction effects become less sizable. Note that although 
the value of mp is below 114 GeV for some points, these 
points can still satisfy the LEP constraints encoded in 
the NMSSMTools because the LEP bound on the MSSM 
Mnp is 92 GeV (see Fig.2 in the paper by Barger et. al. in 
Ref. {2]). By the way, some comprehensive studies on the 
phenomenology of the NMSSM Higgs sector (which may 
be quite different from the phenomenology of the MSSM 
Higgs sector) have been performed in Ref. {2}, where all 
the Higgs couplings and all the decay modes as well as 
productions have been intensively studied. 

From Table[I] we see that although A is scanned in the 
range of —0.5 ~ 0.5, only small values of A (~ 0.1 in 
magnitude) survived the LEP constraints encoded in the 
NMSSMTools. As A gets large, the mixing between sin- 
glet and doublet Higgs fields becomes large (although the 
mixing is not totally determined by À and other param- 
eters are also contributing) and, consequently, is more 
constrained by LEP experiments |2]. Such a small A, to- 
gether with other NMSSM parameters appearing in the 
mass matrix of the Higgs fields, leads to a small mixing 
between singlet and doublet Higgs fields. As shown in 
Table [I] the singlet component UR in A is small, which 
means that the lightest Higgs boson A is not the singlet- 
like one and instead it is doublet-dominant (MSSM-like). 
For the hbb coupling with such a MSSM-like h, the 
NMSSM can still allow for quite different SUSY QCD 
effects compared with the MSSM predictions. The rea- 
son is that, as discussed above, the SUSY QCD effects 
are sensitive to the ratio UH /U#, which is different in 


these two models (even for a small A, the NMSSM can 
still allow for different-size components UŽ and U# in h 
because other NMSSM parameters are also contributing 
to the mass matrix and mixings of the Higgs fields). 

When we switch on the comprehensive experimental 
constraints encoded in the NMSSMTools (10), we obtain 
the results of Asgcp displayed in Fig/4] These experi- 
mental constraints are comprehensive, including the LEP 
II searches for the Higgs boson, various B-decays and the 
muon anomalous magnetic moment (muon g — 2). 
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FIG. 4: Same as Fig. but for ma = 100 GeV in the 
NMSSM. The dark-shaded region denoted by ’x’ (red in 
color) are allowed by the experimental constraints considered 
in the NMSSMTools {10}. The light-shaded region denoted 
by ’e’ (yellow in color) are excluded by LEP experiment or B 
physics, and the medium-shaded region denoted by ’+’ (sky- 
blue in color) are excluded by the muon g — 2. 


From Fig] we see that in the special scenario under 
our consideration (all soft breaking mass parameters are 
degenerate), Msysy is constrained in a certain range. 
The range of Msysy S 600 GeV is not allowed by LEP 
experiments, while Msysy & 1 TeV cannot explain the 
muon g — 2 data (we require the supersymmetric effects 
to account for the deviation of a¢*? — a8/@ = (29.5 + 
8.8) x 10~1° at 2o level). In the allowed region of the 
parameter space, the SUSY QCD corrections to the hbb 
coupling can still be significant, over 20% in magnitude. 
If we switch off the muon g — 2 constraint ( the hadronic 
contribution to a is not so certain aA), the allowed 
parameter space gets much broader, as shown in Fig[4] In 
our study we did not require supersymmetry to explain 
various plausible evidences of dark matter. 

Note that the SUSY contributions to the muon g—2 are 
sensitive to the soft masses in slepton sector and chargino 
sector, but not dependent on the squark or gluino mass 
involved in our SUSY QCD loops. The stringent con- 


straint from the muon g — 2 data shown in Fig[4] comes 
from our simple assumption that all soft masses (in the 
squark sector, the slepton sector and the gaugino sector) 
are degenerate. Of course, the muon g — 2 constraint is 
not necessary; without such an assumption of degeneracy, 
the constraint is lifted. 

Since the SUSY QCD residual effects can be over 20% 
for the hbb coupling and thus over 40% for the ratio of 
the branching fractions (SUSY QCD does not contribute 
to h + 777 at one-loop level) 


Br(h — bb) 


ies ZA T U 
a Br(h— T+T)’ 


(25) 


future measurements at the LHC or ILC may unravel 
such supersymmetric effects. This ratio R,/, is proposed 
in as a probe to new physics. To measure this ratio 
at the LHC, one may count the event numbers of the pro- 
duction of hbb followed respectively by the decay h — bb 
and h — r*7~, and the ratio of these two event-numbers 
can be a measure of R,/, (the difference of efficiency for 
b-tagging and r-tagging should be taken into account). 

Finally, we make some remarks regarding to our re- 
sults: 


(1) We only investigated the SUSY QCD corrections, 
which is (as) and thus should be the most im- 
portant among the SUSY corrections. Among 
the SUSY electroweak corrections, the Higgs-top 


Yukawa corrections, which are O(a), may also 
w 


be sizable although they are seemingly not as large 
as the O(a,) SUSY QCD corrections. With the 
consideration of the SUSY electroweak corrections, 
our conclusion will remain qualatively unchanged. 


(2) Since our main interest is the residual effects 
of heavy sparticles, we assumed all soft SUSY- 
breaking mass parameters are equal to Msysy. 
This is a very strong simplifying assumption. If 
we lift such an assumption and consider the mul- 
tiple free soft parameters, we will obtain the re- 
sults which numerically are different to some ex- 
tent while qualitatively exhibit the same feature, 
i.e., the NMSSM allows broader corrections than 
the MSSM, because the NMSSM parameters can 
complicate the mass matrix and mixings of Higgs 
fields. 


(3) We should stress again that supersymmetry is a 
decoupling theory and the large residual quantum 
effects of sparticles in the Higgs Yukawa couplings 
are present only in case of a light ma. If both 
Msusy and ma take their heavy limit, the SUSY 
effects will vanish. Such a decoupling behavior is 
similar in the NMSSM and MSSM, as shown in 
Fig.2. 


IV. SUMMARY 


We focused on the SUSY QCD effects in the hbb cou- 
pling and performed a comparative study for the two 
popular SUSY models: the MSSM and NMSSM. We 
found that for both models the SUSY QCD can leave 
over large residual quantum effects in hbb coupling if 
the sparticles are relatively heavy (a few TeV) while the 
Higgs sector is not so heavy (m4 is not so large). Com- 
pared with the MSSM results, the NMSSM can allow 
for a much broader region of such residual effects. Since 
these residual effects can be over 20% in magnitude, fu- 
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